Comment on "On Novel attractive forces between ions in quantum plasmas — failure 

of linearized quantum hydrodynamics" 
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In a recent paper Bonitz, Pehlke and Schoof hereafter referred to as BPS, have raised some 
points against the newly found Shukla-Eliasson attractive potential @, @] , hereafter refererred to as 
SEAP, around a stationary test charge in a quantum plasma Our objective here is to discuss 

the inappropriateness of BPS reasoning concerning the applicability of the linearized quantum hy- 
drodynamic theory, as well as to point out the shortcomings in BSP's arguments and to suggest for 
rescuing the BSP's density functional theory and simulations which have failed to produce results 
that correctly match with that of Shukla and Eliasson @, 0|- 
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Bonitz, Pehlke and Schoof 'l| (BPS) have compared 
the findings from their density functional theory (DFT) 
simulations with that of an analytical theory of Shukla- 
Eliasson (SE) dm which reported a short-range (at the 
scale size of several Bohr atomic radii) attractive poten- 
tial (AP) distribution around a stationary ion test charge 
in a quantum plasma. To start with, we need to evalu- 
ate the assumptions in the quantum plasma regime em- 
ployed in the theory of Shukla and Eliasson (SE). SE 
Q have used the well established 0, [f| quantum hy- 
drodynamic (QHD) model for degenerate electrons to 
calculate the electric potential around a stationary ion 
test charge by using the electron dielectric constant de- 
duced from the linearized electron continuity nonrela- 
tivistic electron momentum and Poisson equations. The 
electron momentum equation includes the electrostatic 
and quantum forces (e.g the non-relativistic quantum 
statistical electron pressure in the zero-temperature free- 
electron Fermi gas limit, electron-exchange and electron 
correlation effects due to the electron spin, and the quan- 
tum electron recoil effect associated with overlapping of 
electron wave functions and their dispersion at atomic 
scales). Thus, the QHD model is subject to some ap- 
proximations regarding quantum degeneracy of free elec- 
trons in a Fermi plasma and comparison of the quantum 
statistical pressure and the quantum recoil effect associ- 
ated with electron tunneling through the quantum Bohm 
potential. The quantum electron degeneracy comes into 
the picture when the thermal de Broglie wavelength of 
electrons, Xb — h/ ^kBT p m e is comparable with the av- 
erage inter-electron spacing tq k, (3/47rno) 1 ^ 3 , in addi- 
tion to the condition Xb <C Xl, where h is the Planck 
constant divided by 2ir, ks the Boltzmann constant, 
T p the plasma temperature, m e the rest mass of elec- 
trons, no the electron number density, = e 2 /ksT p 



the Landau length [6|, and e the magnitude of electron 
charge. However, the free electron assumption, which is 
implicitly hidden in the Fermi-Dirac statistical degener- 
acy pressure Pdeg ,7], is often overlooked in literature. 
It is important to note that the latter assumption is 
only valid for metallic compounds with very large elec- 
trical conductivity 0. This condition is also valid for 
solid density quantum plasmas and warm dense matter 
(WDM). Thus, the SEAP theory critically relies on the 
free-electron assumption, in addition to the zero-Fermi 
temperature electron degeneracy with (huj p /kBTp) 2 ~ 1, 
and any application of this theory beyond those assump- 
tions can lead to nonphysical results, as discussed be- 
low. Here lo p — (iirnoe 2 /rrie) 1 ^ 2 is the electron plasma 
frequency and Tp = (h 2 /2m e fcs)(37r 2 no) 2 ^ 3 the Fermi 
electron temperature. The origin of the negative attrac- 
tive SE potential is attributed to the consideration of 
the quantum recoil effect, which may dominate over the 
quantum statistical electron pressure as well as electron- 
exchange and electron-correlation effects, depending on 
the certain plasma density ranges discussed in Ref. [3j. 

BPS have compared the SEAP theory with their DFT 
simulations for a hydrogen plasma. Their code reveals 
interesting mismatch between the two theories and there- 
fore they concluded that the linearized QHD theory fails. 
Although the DFT theory is far from exact 0| and still 
subject to many optimizations (Tol j with many different 
version of it found in the literature, one should seek the 
cause of fundamental divergences between the two ap- 
proaches (the DFT and QHD models) in the inadequate 
comparison. BPS throughout their PRE paper [l| fre- 
quently refer to Friedel oscillations as a possible mech- 
anism for the SE effect. We refute to such an assertion 
because, Friedel oscillations arise from localized pertur- 
bations in metallic and semiconducting materials and are 
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absent in insulators with molecular bindings. Friedel os- 
cillations are closely related to the Kohn anomaly 11 1 



or even r, 



caused by singularity in the plasma wave dispersion rela- 
tion [l2| (regardless of the quantum recoil effect at atomic 
scales) at the wavenumber value of k = 2kp, where fcp 
is the Fermi wavenumber of electrons. Lindhard and Vi- 
densk [l2[ have obtained the dielectric constant for the 
free electron liquid in metals, which becomes singular at 
some k values. Although, such singularity is negligible 
in reciprocal space, however, its Fourier transform or its 
image in the real space causes a strong oscillations in the 
real-part of the dielectric function in proximity of the 
singularity, due to the well known Gibbs phenomenon. 
Such oscillations are reflected in the screenin g p otential 
of the test charge in a Fermi electron liquid [13(. How- 
ever, no such oscillations (Friedel oscillation of same ori- 
gin) has been reported for quantum plasmas in which 
the Fermi- wavevector and the Fermi surface are not well- 
defined parameters. In quantum plasmas without the 
quantum recoil effect, one encounters only the Thomas- 
Fermi (TF) short-range [of the order of the TF screening 
radius (ksTp jm e ) x l 2 /oj p ] repulsive potential. 

Now, let us closely inspect the points raised by BPS 
against the SEAP. In their Fig. 1, BPS [1] have shaded 
a region, they call the region of validity of the linearized 
QHD, for which hujn < ksTp. However, the QHD model 
is actually valid [3] for (hui p /k B T F ) 2 ~ 1. There- 
fore, the shaded area in Fig. 1 of BPS does not apply 
to dense plasmas with zero-temperature completely de- 
generate electron fluids where electron-ion collisions are 
greatly inhibited by the Pauli blocking mechanism. This 
condition may be relaxed to the quantum coupling pa- 
rameter, gQ — (hLUp/kBTp) 2 [ , even slightly greater than 
unity in the weak coupling limit [5]. In fact, the weak 
coupling limit in quantum plasmas is not a strict one and 
it closely relates to the mean-free-path of electrons which 
may be very large in the quantum plasma cases like some 
good conductors. In particular, it is well-known that in 
metallic compounds the characteristic plasmon frequency 
is at least one order of magnitude larger than that of 
electron-ion and many orders of magnitudes larger than 
that for electron-electron collisions [5|]. As it will be ap- 
parent bellow, it is however irrelevant to apply the QHD 
degenerate free electron model used in SEAP theory 
to nonmetallic hydrogen densities, characterized by the 
Brueckner parameter (r s = tq/tb with tb — h 2 /e 2 m e 
being the Bohr radius), lying in the region gQ > 1. 

BPS in their Fig. 5 [l[ compare their DFT simulations 
with SEAP for three different Brueckner parameters of 
r s = 7,4,1.5 corresponding to electron number densi- 
ties of n ~ (4.7 x 10 21 ,2.5 x 10 22 ,4.8 x 10 23 ) citT 3 , 
respectively. Since, the existence of the SEAP is a di- 
rect consequence of the interplay between the quantum 
forces which strongly rely on the free electron model and 
complete degeneracy assumptions, it is readily apparent 
that the SEAP does not apply for hydrogen with r s = 7 



4. In fact, BPS correctly remark in p. 3 
of their paper [l[ that in the range r s > 3 hydrogen gas 
is in its nonmetallic molecular binding state and these 
bounds break only for r s ~ 2 ... 3. That is, what they 
compare in their Fig. 5 for r s — 4 and r s — 7 is irrel- 
evant for the SEAP, since it overrides the required as- 
sumptions for the applicability of the linearized QHD. 
In other words, the SEAP theory is an strongly ionized 
atomic theory applied only to highly conductive (pres- 
sure) ionized materials with densities beyond the Mott 
metal-insulator transition, which happens to be in the 
regime r s = 1.2... 1.5 (e.g. see p. 3 of Ref. [l[). In 
fact, it is clearly observed from Fig. 5 of BPS that the 
best match between the two theories coincides with the 
value of r s — 1.5, which is close to the Mott transition 
for hydrogen composition. Therefore, it should be kept 
in mind that, despite the apparent wide density range for 
which the SEAP minimum exists in hydrogen plasma, the 
SEAP unlike the DFT is neither a theory for description 
of molecular bindings nor it gives rise to the molecular 
Lennard- Jones potential. Therefore, one must critically 
examine the applicability of the free electron assumption, 
which can only be valid for a metallic pressure-ionized 
hydrogen in the limit r s < 1.2. Obviously, this does not 
count as a defect for the SEAP theory, but, it is left 
for the reader to examine the outcome against the basic 
assumptions for the QHD model. 

The Shukla-Eliasson attractive potential @ , which can 
lead ion-ion correlations, is caused by electron density lo- 
calization due to quantum electron dispersion/recoil ef- 
fect. It can be observed from Fig. 5 of BPS that for den- 
sities of r s = 1 . . . 0.6 the SEAP minimum can become 
more pronounced and the minimum of BPS disappears. 
Such defect in the standard DFT simulations may be 
overcome by extending the conventional Thomas-Fermi 
(TF) screening of ions to the modified Thomas-Fermi- 
Weizsacker (TFW) improved model, which also includes 
the quantum electron recoil effect [10( . It has been shown 
(e.g. see Fig. 6 of Ref. [ijj) that the improved TFW 
model can give rise to much deeper potential valley than 
that calculated with the ordinary TF model. The DFT 
theory with improved TFW density model has been also 
found to be more consistent with laboratory data com- 
pared to the standard DFT theories. Furthermore, BPS 
in their paper call the DFT simulation method as a ref- 
erence to other theoretical results. This is by the way far 
from exact, since, there are major complexities associated 
with DFT approximate calculations [15[ , which should be 
overcome, before it can be claimed as a reference theory. 

In summary, the material presented in BPS paper is 
partially (where related to molecular binding) irrelevant 
to the SEAP theory. On the other hand, BPS in their 
PRE paper [l[ have erroneously attributed the existence 
of the SEAP minimum to the Friedel oscillations (related 
to the singularity in the Lindhard dielectric constant for 
ordered Fermi liquids possessing sharp Fermi surface), 
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relevant to metals and semiconductors with a well denned 
Fermi- wavevector. The SEAP minimum has a direct root 
in the electron quantum recoil effect and electron wave- 
function interferences, present in any quantum plasma in 
the absence of a well-defined Brillouin zone. Moreover, 
the free-electron model assumption used in the SEAP 
theory has been overlooked by BPS in their PRE paper. 
The DFT theory, in its current state, is far from being 
complete in the context of plasma physics with an en- 
semble of degenerate electrons interacting in a collective 
fashion at atomic dimensions, and many improvements 
are under way and until then one is unable to claim DFT 
as an ultimate theory. Finally, the extent of validity of 
the linear hydrodynamic description of ion structure fac- 
tors in dense plasmas has been investigated by Mithen 
et al. fl6j ] who concluded that such an approach can be 
used to effectively model the ion response in compressed 
plasmas for a wide range of the plasma number densities 
that can be probed experimentally. 
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